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We study the phenomenon of length scale competition, an instability of solitons and other coher- 
ent structures that takes place when their size is of the same order of some characteristic scale of 
the system in which they propagate. Working on the framework of nonlinear Klein-Gordon models 
as a paradigmatic example, we show that this instability can be understood by means of a collective 
coordinate approach in terms of soliton position and width. As a consequence, we provide a quanti- 
tative, natural explanation of the phenomenon in much simpler terms than any previous treatment 
of the problem. Our technique allows to study the existence of length scale competition in most 
soliton bearing nonlinear models and can be extended to coherent structures with more degrees of 
freedom, such as breathers. 
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LEAD PARAGRAPH 

Solitons, solitary waves, vortices and other co- 
herent structures possess, generally speaking, a 
characteristic length or size. One important fea- 
ture of these coherent structures, which usually 
are exact solutions of certain nonlinear models, 
is their robustness when the corresponding mod- 
els are perturbed in different ways. A case rel- 
evant in many real applications is that of space- 
dependent perturbations, that may or may not 
have their own typical length scale. Interestingly, 
in the latter case, it has been known for a decade 
that when the perturbation length scale is com- 
parable to the size of the coherent structures, the 
effects of even very small perturbing terms are 
dramatically enhanced. Although some analyti- 
cal approaches have shed some light of the mecha- 
nisms for this special instability, a clear-cut, sim- 
ple explanation was lacking. In this paper, we 
show how such explanation arises by means of a 
reduction of degrees of freedom through the so- 
called collective coordinate technique. The an- 
alytical results have a straightforward physical 
interpretation in terms of a resonant-like phe- 
nomenon. Notwithstanding the fact that we work 
on a specific class of soliton-bearing equations, 
our approach is readily generalizable to other 
equations and/or types of coherent structures. 
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I. INTRODUCTION 

Fifty years after the pioneering discoveries of Fermi, 
Pasta and Ulam Q , the paradigm of coherent structures 
has proven itself one of the most fruitful ones of Nonlinear 
Science (2i- Fronts, solitons, solitary waves, breathers, 
or vortices are instances of such coherent structures of 
relevance in a plethora of applications in very different 
fields. One of the chief reasons that gives all these nonlin- 
ear excitations their paradigmatic character is their ro- 
bustness and stability: Generally speaking, when systems 
supporting these structures are perturbed, the structures 
continue to exist, albeit with modifications in their pa- 
rameters or small changes in shape (see for reviews). 
This property that all these objects (approximately) re- 
tain their identity allows one to rely on them to interpret 
the effects of perturbations on general solutions of the 
corresponding models. 

Among the different types of coherent structures one 
can encounter, topological solitons are particularly ro- 
bust due to the existence of a conserved quantity named 
topological charge. Objects in this class are, e.g., kinks 
or vortices and can be found in systems ranging from 
Josephson superconducting devices to fluid dynamics. A 
particularly important representative of models support- 
ing topological solitons is the family of nonlinear Klein- 
Gordon equations whose expression is 

0tt - 4>xx + ^ = 0- (1) 
d(j) 

Specially important cases of this equation occur when 
[/((/)) = j{(tP — 1)^, giving the so-called (j)'^ equation, and 
when U{(t)) = 1 — cos(0), leading to the sine-Gordon (sG) 
equation, which is one of the few examples of fully in- 
tegrable systems 0- Indeed, for any initial data the 



2 



solution of the sine-Gordon equation can be expressed 
as a sum of kinks (and antikinks), breathers, and linear 
waves. Here we focus on kink solitons, which have the 
form 



/ X — vt 



(f>{x, t) = 4 arctan < exp . 

I V V 1 



(2) 



< V < 1 being a free parameter that specifies the kink 
velocity. The topological character of these solutions 
arises from the fact that they join two minima of the 
potential U {(/)), and therefore they cannot be destroyed 
in an infinite system. Our other example, the 0** equa- 
tion, is not integrable, but supports topological, kink-like 
solutions as well, given by 



(x, t) = tanh 



X — vt 



(3) 



It is by now well established, already from pioneering 
works in the seventies 6, 7] that both types of kinks be- 
have, under a wide class of perturbations, like relativis- 
tic particles. The relativistic character arises from the 
Lorentz invariance of their dynamics, see Eq. and 
implies that there is a maximum propagation velocity 
for kinks (1 in our units) and their characteristic width 
decreases with velocity. Indeed, even behaving as parti- 
cles, kinks do have a characteristic width; however, for 
most perturbations, that is not a relevant parameter and 
one can consider kinks as point-like particles. This is not 
the case when the perturbation itself gives rise to certain 
length scale of its own, a situation that leads to the phe- 
nomenon of length-scale competition, first reported in Q 
(see for a review). This phenomenon is nothing but 
an instability that occurs when the length of a coherent 
structure approximately matches that of the perturba- 
tion: Then, small values of the perturbation amplitude 
are enough to cause large modifications or even destroy 
the structure. Thus, in |8j, the perturbation considered 
was sinusoidal, of the form 

(jjtt- (t)xx + ^{l + ecos{kx)) = 0, (4) 

where e and k are arbitrary parameters. The structures 
studied here were breathers, which are exact solutions of 
the sine-Gordon equation with a time dependent, oscilla- 
tory mode (hence the name 'breather') and that can be 
seen as a bound kink-antikink pair. It was found that 
small k values, i.e., long perturbation wavelengths, in- 
duced breathers to move as particles in the sinusoidal po- 
tential, whereas large k or equivalent short perturbation 
wavelengths, were unnoticed by the breathers. In the in- 
termediate regime, where length scales were comparable, 
breathers (which are non topological) were destroyed. 

As breathers are quite complicated objects, the issue 
of length scale competition was addressed for kinks in 
[T(tI |. In this case, kinks were not destroyed because of 
the conservation of the topological charge, but length 



scale competition was present in a different way: Keep- 
ing all other parameters of the equation constant, it was 
observed that kinks could not propagate when the per- 
turbation wavelength was of the order of their width. In 
all other (smaller or larger) perturbations, propagation 
was possible and, once again, easily understood in terms 
of an effective point-like particle. Although a explana- 
tion of this phenomenon was provided in in terms 
of a (numerical) linear stability analysis and the radia- 
tion emitted by the kink, it was not a fully satisfactory 
argument for two reasons: First, the role of the kink 
width was not at all transparent, and second, there were 
no simple analytical results. These are important issues 
because length scale competition is a rather general phe- 
nomenon: It has been observed in different models (such 
as the nonlinear Schrodinger equation or with other 
perturbations, including random ones |l2l |. Therefore, 
having a simple, clear explanation of length scale compe- 
tition will be immediately of use in those other contexts. 

The aim of the present paper is to show that length 
scale competition can be understood through a collec- 
tive coordinate approximation. Collective coordinate ap- 
proaches were introduced in |^ Q to describe kinks as 
particles (see ^^AiM a very large number of differ- 
ent techniques and applications of this idea). Although 
the original approximation was to reduce the equation 
of motion for the kink to an ordinary differential equa- 
tion for a time dependent, collective coordinate which 
was identified with its center, it is being realized lately 
that other collective coordinates can be used instead of 
or in addition to the kink center. One of the most natural 
additional coordinates to consider is the kink width, an 
approach that has already produced new and unexpected 
results such as the existence of anomalous resonances 
or the rectification of ac drivings There are also 

cases in which one has to consider three or more collec- 
tive coordinates (see, e.g., ^3)- only natural then 
to apply these extended collective coordinate approxima- 
tions to the problem of length scale competition, in search 
for the analytical explanation needed. As we will see be- 
low, taking into account the kink width dependence on 
time is indeed enough to reproduce the phenomenology 
observed in the numerical simulations. Our approach is 
detailed in the next Section, whereas in Sec. 3 we collect 
our results and discuss our conclusions. 



II. COLLECTIVE COORDINATE APPROACH 

We now present our collective coordinate approach to 
the problem of length scale competition for kinks [lol| . 
We will use the lagrangian based approach developed in 
|16| . which is very simple and direct. Equivalent results 
can be obtained with the so-called generalized travelling 
wave Ansatz 17], somewhat more involved in terms of 
computation but valid even for systems that cannot be 
described in terms of a lagrangian. 

Let us consider the generically perturbed Klein- 
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Gordon equation: 



a) 



0tt-0x:r + ^+e/(a;,i)5(0) = 0. (5) 

The starting point of our approach is the lagrangian for 
the above equation, which is given by 



200 r 



<^9{(l>)4>x / dyf{y,t) 



(6) 



As stated above, we now focus on the behaviour of kink 
excitations of the form Q and To do so, we wiU use 
a two collective coordinates approach by substituting the 
Ansatz 

(t){z{t)) =teinh{z{t)) (7) 
in the lagrangian of the 0^ system, and 

0(z(t)) = 4arctan{exp(z(t))} (8) 

in sG, where z{t) — and X{t) and l{t) are two 

collective coordinates that represent the position of the 
center and the width of the kink, respectively. Substitut- 
ing the expresions Q and in the lagrangian ® with 



dU 
d4> ■ 



we 



our perturbation, f{x,t) = cos{kx) and g{4>) 
obtain an expresion of the lagrangian in terms of X and 
I, 



MqIo 
21 



aMoloj2 _Ml(b.j^L 
21 2 \l Iq 



-- cos{kX)w{a)\a=ki, 



where, for the (f)"^ system, AIq — 4/(3\/2), = 
a = (tt^ — 6)/12, and for the sG system, Afo = 

^2 , 



and a — -k jVl, and 



dz tanh(z)(0'(z))^ sin(az) 



which is 



for (h and 



w{a) 



-KO^ia^ + 4) 
24sinh(^) 



(9) 
V2 and 



(10) 
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FIG. 1: Different behaviours of X(i) for different values of 
fc for e = 0.7, X(0) = 0, X(0) = 0.5, /(O) = /o(l - XiS^ffl'^ 
and Z(0) =0 in the a) 0* model, b) sG model. 



The equations of motion of X and / can now be ob- 
tained using the Lagrange equations, 



d_ [dL_\ _ dL_ 



(13) 



where Yi stands for the collective coordinates X and I. 
The ODE system for X and I is, finally. 



P = esmikX)w{a)\a=kh 
a 



Q = 



MM ( 1 



2Mo«o 

tco%{kX)w' {a)\a=ki, 



12 12 

L Iq 



(14) 



(15) 



w{a) 



2na^ 



(12) 



for the sG model. We note that the effect of a spatially 
periodic perturbation like the one considered here was 
studied for the (f)'^ model in although the authors 
were not aware of the existence of length scale competi- 
tion in this system and focused on unrelated issues. 



where P = MqIoX/I and Q = aMM/l. 

The equations above are our final result for the dy- 
namics of sG and (j)^ kinks in terms of their center and 
width. As can be seen, they are quite complicated equa- 
tions and we have not been able to solve them analyti- 
cally. Therefore, in order to check whether or not they 
predict the appearance of length scale competition, we 
have integrated them numerically using a Runge-Kutta 
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FIG. 2: Different behaviours of for different values of k 
for e = 0.7, X(0) = 0, X(0) = 0.5 as obtained from numerical 
simulations of the full 0* model. 



scheme [19(. The simulation results for different values 
of k, X(0) = are in Fig. [T] We have taken e = 0.7 
in order to conipare with the numerical results for the 
sG equation in [13 . The plots present already the phys- 
ical variable i.e., the position of the center of the 
kink. As we may see, the agreement with the results in 
[lOl is excellent for the sG equation (cf. Figs. 7 and 8 
in lOj): For small and large wavelengths the kink, ini- 
tially at the top of one of the maxima of the perturba- 
tion, travels freely; however, for intermediate scales, the 
kink is trapped and even moves backward for a while. As 
previous works 18] on the same problem for the 0"* equa- 
tion did not touch upon this issue, we have carried out 
numerical simulations, also using a Runge-Kutta scheme, 
of the full partial differential equation. The results, pre- 
sented in Fig. [51 confirm once again the agreement with 
the prediction of the collective coordinate approach. We 
see that for /c = 2, the value at which the collective coor- 
dinate prediction is more dramatic, the kink is trapped 
at the first potential well, indicating that for that value 
the length scale competition phenomenon is close to its 
maximum effect. Interestingly, a more detailed analysis 
of Fig. 121 shows two different kink trapping processes. In 
agreement with the collective coordinate prediction, for 
/c = 2, 3 we see that the kink is trapped very early, af- 
ter having travelled a few potential wells at most. This 
is the length scale competition regime. However, there 
is an additional trapping mechanism: As is well known 
|18|. kinks travelling on a periodic potential emit radi- 
ation. This process leads to a gradual slowing down of 
the kink until it is finally unable to proceed over a fur- 
ther potential well. This takes place at a much slower 
rate than length scale competition and hence the kink is 
stopped after travelling a larger distance in the system. 



The radiation emission as observed in the simulation is 
smoother than in the length scale competition trapping. 



III. DISCUSSION AND CONCLUSIONS 

As we have seen in the preceding section, a collec- 
tive coordinate approach in terms of the kink center and 
width is able to explain in a correct quantitative manner 
the phenomenon of length scale competition, observed in 
numerical simulations earlier for the sG equation with a 
spatially periodic perturbation |^ ^| . The structure of 
the equations makes it clear the necessity for a second 
collective coordinate; imposing l(t) — Iq constant, we re- 
cover the equation for the center already derived in [lol |. 
which shows no sign at all of length scale competition, 
predicting effective particle-like behavior for all k. The 
validity of this approach has been also shown in the con- 
text of the (f)'^ equation, which had not been considered 
before from this viewpoint. In spite of the fact that the 
collective coordinate equations cannot be solved analyti- 
cally, they provide us with the physical explanation of the 
phenomenon in so far as they reveal the key role played 
by the width changes with time and their coupling with 
the translational degree of freedom. 

It is interesting to reconsider the analysis carried out 
in [10| of length scale competition through a numerical 
linear stability analysis. In that previous work, it was ar- 
gued that the instability arose because, for the relevant 
perturbation wavelengths, radiation modes moved below 
the lowest phonon band, inducing the emission of long 
wavelength radiation which in turn led to the trapping 
of the kink. It was also argued that those modes became 
internal modes, i.e., kink shape deformation modes in the 
process. The approach presented here is a much more 
simple way to account for these phonon effects: Indeed, 
as was shown by Quintero and Kevrekidis |23|, (odd) 
phonons do give rise to width oscillations very similar to 
those induced by an internal mode [23 . We are confident 
that what our perturbation technique is making clear is 
precisely the result of the action of those phonons, sum- 
marized in our approach in the width variable l{t). The 
case for the 0^ equation is sligthly different: Whereas 
the sG kink does not have an internal mode and l{t) 
is hence a collective description of phonon modes, the 0^ 
kink possesses an intrinsic internal mode that is easily ex- 
cited by different mechanisms (such as interaction with 
inhomogeneities, [l^,!!^)- Therefore, one would expect 
that for the cj)'^ equation the effect of a perturbation of a 
given length is more dramatic than for the sG model, as 
indeed is the case: See Fig. ^ for a comparison, showing 
that the kink is trapped for a wider range of values of k 
in the 0"* equation. Remarkably, the present and previ- 
ous results using this two collective coordinate approach, 
and particularly their interpretation in terms of phonons, 
suggest that this technique could be something like a 'sec- 
ond order collective coordinate perturbation theory', the 
width degree of freedom playing the role of second order 
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term. We believe it is appealing to explore this possibil- 
ity from a more formal viewpoint; if this idea is correct, 
then one could think of a scheme for adding in an stan- 
dard way as many collective coordinates as needed to 
achieve the required accuracy. Progress in that direction 
would provide the necessary mathematical grounds for 
this fruitful approximate technique. 

Finally, some comments are in order regarding the ap- 
plicability of our results. We believe that the collective 
coordinate approach may also explain the length scale 
competition for breathers, that so far lacks any expla- 
nation. For this problem, the approach would likely in- 
volve the breather center and its frequency, as this mag- 
nitude controls the kink width when the kink is at rest. 
Of course, such an Ansatz would only be valid for the 
breather initially at rest, and the description of the dy- 
namical problem would be more involved, needing per- 
haps more collective coordinates (such as an independent 
width). If this approach succeeds, one can then extend 
it to other breather like excitations, such as nonlinear 
Schrodinger solitons or intrinsic localized modes. Work 
along these lines is in progress. On the other hand, it 



would be very important to have an experimental setup 
where all these conclusions could be tested in the real 
world. A modified Josephson junction device has been 
proposed recently '23| where the role of the kink width 
is crucial in determining the performance and dynam- 
ical characteristics. We believe that a straightforward 
modification of that design would permit an experimen- 
tal verification of our results. We hope that this paper 
stimulates research in that direction. 
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